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Abstract 



I give a pedagogical summary of the methods of light-cone variables, rapidity 
and pseudo-rapidity. Then I show how these methods are useful in analyzing 
diffractive high-energy collisions. 



I. INTRODUCTION 



The use of light-cone variables, rapidity and pseudo-rapidity is very common in treating 
high-energy scattering, particularly in hadron-hadron and lepton-hadron collisions. The 
essential features of these collisions that make these variables of utility are the presence of 
ultra-relativistic particles and a preferred axis. 

In these notes I will explain the theory of these variables. The only prerequisites are 
a knowledge of the elements of special relativity and an acquaintance with the general 
phenomenological features of high-energy scattering. Then I will show how these variables 
can be used conveniently in the analysis of the kinematics of diffractive processes. This will 
include the derivation of an approximate relation between the size of the rapidity gap and 
the variable Xf. 

The material in these notes is essentially all standard and well known to many workers 
in the field. But it is not always easy to find the material conveniently presented. 



Light-cone coordinates are defined by a change of variables from the usual (t, x, y, z) (or 
(0, 1,2,3)) coordinates. Given a vector V^, its light-cone components are defined by 



II. LIGHT-CONE COORDINATES, RAPIDITY 



A. Definitions 



V 
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and I will write = {V~^,V , V-^). Some authors prefer to omit the 1/a/2 factor in Eq. 
(0). It can easily be verified that Lorentz invariant scalar products have the formQ 

V-W = V+W- + V'W+ - Vt ■ Wt, 

V-V = - VS- (2) 

What are the motivations for defining such coordinates, which evidently depend on a 
particular choice of the z axis? One is that these coordinates transform very simply under 
boosts along the 2;-axis. Another is that when a vector is highly boosted along the z axis, 
light-cone coordinates nicely show what are the large and small components of momentum. 
Typically one uses light-cone coordinates in a situation like high-energy hadron scattering. 
In that situation, there is a natural choice of an axis, the collision axis, and one frequently 
needs to transform between different frames related by boosts along the axis. Commonly 
used frames include the rest frame of one of the incoming particles, the overall center-of-mass 
frame, and the center-of-mass of a partonic subprocess. 



B. Boosts 

Let us boost the coordinates in the z direction to make a new vector V'^. In the ordinary 
{t, X, y, z) components we have the well known formulae 

T/O _|_ .,yz „,T/0 _|_ yz 

y>o ^ V y,^ ^ ^ v"" = v'y = yy. (3) 



V 



It is easy to derive the following for the light-cone components: 

= yi- = y-Q-^^ y't = Vt, (4) 

where the hyperbolic angle is | In so that v = tanh-?/'. 

Notice that if we apply two boosts of parameters ipi and 1^2 the result is a boost tpi+ijj2- 
This is clearly simpler than the corresponding result expressed in terms of velocities. 



C. Rapidity 

a. Boost of particle momentum Consider a particle of mass m that is obtained by a 
boost from the rest frame. Its momentum is 

2 

m 



I prefer to make a distinction between contravariant vectors, whose indices are superscripts, and 
covariant vectors, whose indices are subscripts. In all contexts where the distinction matters I will 
put the Lorentz indices in their correct places. 
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Notice that if the boost is very large (positive or negative), only one of the two non-zero 
light-cone components of p'^ is large; the other component becomes small. With the usual 
components two of the components, p° and p^, become large. 

Suppose next that we have two such particles, pi and p2, with the boost for particle 1 
being much larger than that for particle 2. Then in the scalar product of the two momenta 
only one component of each momentum dominates the result, thus for example {P1+P2Y — 
2piP2 ■ This implies that, when analyzing the sizes of scalar products of highly boosted 
particles, it is simpler to use light-cone components than to use conventional components. 

b. Definition of rapidity Since the ratio p'^ /p^ gives a measure e^'^ of the boost from 
the rest frame, we are led to the following definition of a quantity called "rapidity" : 

1, p+ 1, E + p' 
2 p 2 h — p^ 

which can be applied to a particle of non-zero transverse momentum. The 4-momentum of 
a particle of rapidity y and transverse momentum is 



^=|./!!l!±iie»,,/!!l!±»e-».pH, (7) 



with a/ m? + p\ being called the transverse energy Et of the particle. It can be checked 
that the scalar product of two momenta is 

pi • P2 = E1TE2T cosh(yi - ^2) - Pit • P2r- (8) 

In the case where the transverse momenta are negligible, this reduces to m^cosh(yi — 1/2), 
which is like the formula for the product of two Euclidean vectors pi • P2 = P1P2 cos 6, with 
the trigonometric cosine being replaced by the hyperbolic cosine. 

c. Transformation under boosts Under a boost in the z direction, rapidity transforms 
additively: 

y^y'^y + 'il^- (9) 

This implies that in situations where we have a frequent need to work with boosts along the 
z axis it is economical to label the momentum of a particle by its rapidity and transverse 
momentum, rather than to use 3-momentum. 

d. Rapidity distributions in high-energy collisions It also happens that in most collisions 
in high-energy hadronic scattering, the distribution of final-state hadrons is approximately 
uniform in rapidity, within kinematic limits. That is, the distribution of final-state hadrons 
is approximately invariant under boosts in the z direction. This implies that rapidity and 
transverse momentum are appropriate variables for analyzing data and that detector ele- 
ments should be approximately uniformly spaced in rapidity. (What is physically possible 
is to make a detector uniform in the angular variable pseudo-rapidity that I will discuss 
below.) This is in contrast to the situation for e~^e~ collisions where most of the interest 
is in events generated via annihilation into an electro-weak boson. Such events are much 
closer to uniform in solid angle than uniform in rapidity. 
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e. Non-relativistic limit Observe that for a non-relativistic particle, rapidity is the same 
as velocity along the ^-axis, for then 



y = - Xyv — ^ - In = V . (10) 

2 E ~ 2 m — mv^ 

Non-relativistic velocities transform additively under boosts, and the non-linear change of 
variable from velocity to rapidity allows this additive rule Eq. (^ to apply to relativistic 
particles (but only in one direction of boost). 

One way of seeing this is as follows: The relativistic law for addition of velocities in one 
dimension is 

where /3i2 is the velocity of some object 1 measured in the rest-frame of object 2, etc. This 
formula is reminiscent of the following property of hyperbolic tangents: 

, / . tanh A + tanh _B 

tanh(A + B) = . (12) 

^ ^ l + tanhAtanhE ^ ^ 

So to obtain a linear addition law, we should write f3i2 = tanhyli2, and then the rule Eq. 
(|rT|) for the addition of velocities becomes simply Ai^ = A12 + A23. The A variables are 
exactly relative rapidities, since 

v' = ^ = = tanh y. (13) 

f. Relative velocity Rapidity is the natural relativistic velocity variable. Suppose we 
have a proton and a pion with the same rapidity at pt = 0. Then they have no relative 
velocity; to see this, one just boosts to the rest frame of one of the particles. But these same 
particles have very different energies: Ep = ^E^^. 

D. Pseudo-rapidity 

As I will now explain, the rapidity of a particle can easily be measured in a situation 
where its mass is negligible, for then it is simply related to the polar angle of the particle. 
First let us define the pseudo-rapidity of a particle by 

g 

7] = — In tan -, (14) 

where 6 is the angle of the 3-momentum of the particle relative to the +z axis. It is easy to 
derive an expression for rapidity in terms of pseudo-rapidity and transverse momentum: 



'm^ + p'i, cosh^ + Pt sinh 77 

In-^^ ^ . (15) 

^Jm? + p\ 
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In the limit that m ^ px, y ^ rj. This accounts both for the name 'pseudo-rapidity' and 
for the ubiquitous use of pseudo-rapidity in high-transverse-momentum physics. Angles, 
and hence pseudo-rapidity, are easy to measure. But it is really the rapidity that is of 
physical significance: for example the distribution of particles in a minimum bias event is 
approximately uniform in rapidity over the kinematic range available. 

The distinction between rapidity and pseudo-rapidity is very clear when one examines 
the kinematic limits on the two variables. In a collision of a given energy, there is a limit 
to the energy of the particles that can be produced. This can easily be translated to limits 
on the rapidities of the produced particles of a given mass. But there is no limit on the 
pseudo-rapidity, since a particle can be physically produced at zero angle (or at 180°), where 
its pseudo-rapidity is infinite. The particles for which the distinction is very significant are 
those for which the transverse momentum is substantially less than the mass. Note: from 
Eq. (|15D it follows that y < r] always. 



III. DIFFRACTIVE SCATTERING 
A. Kinematics 

Let us define a single diffractive collision as one in which one of the beam particles 
emerges almost unscathed from the scattering, with only a small deflection and a small loss 
of energy (Fig. |T|): 

p + p p' + X. (16) 

For the sake of definiteness, I have assumed a proton-antiproton collision with the proton 
being diffracted. One can treat any other collision in the same fashion. Indeed, replacing 
the incoming and outgoing protons by electrons gives us electroproduction processes, and 
the low events that are used to measure photoproduction cross sections have exactly the 
same kinematic properties as diffractive hadron collisions. 

Suppose the particle p is moving in the +z direction, and p in the opposite direction. 
Then we define the kinematics in terms of the transverse momentum p'^ of the diffracted 
hadron and a variable xp defined as 1 — p'^ jp'^ . Thus the components of p and p' are 

""K'^-^^'-^f?^-'''^)- '''' 

It is convenient to examine the momentum transfer: 

Clearly, xp represents a fractional loss of the -|- component of momentum. When the two 
protons are both moving fast, this fraction is also the fractional energy loss: Xp = 1 — 
E' / E+0{1/ E)] but the definition in light cone components converts this to a boost-invariant 
definition. The invariant momentum transfer is 
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FIG. 1. Diffr active scattering. 



1^ 2 2 

^ _ ^2 _ P T + 

1 - Xp 

—p% as xp 0. (19) 

Thus at small xp, t is an equivalent variable to the transverse momentum of the outgoing 
diffracted hadron. 

The diffractive region is when xp is small enough that the power law attributed to 
pomeron exchange dominates. The exact nature of the pomeron in QCD is still quite 
uncertain. But from a practical point-of-view one can identify it with whatever is exchanged 
between hadrons to make the approximately constant total cross sections that are observed 
at high-energy. The object exchanged between the protons in Fig. |I] and the system labeled 
X is a pomeron when xp is sufficiently small. In the Regge phenomenology of diffractive 
scattering, this is associated with an xp distribution of the events that is approximately 
dxp/xp. 

A more general definition of diffraction is that it consists of events in which there are non- 
exponentially suppressed rapidity gaps. The events that are diffractive in the first definition 



in fact are also diffractive according to the second definition. As I will show in Sect. P^IIE| , 
events with small xp have a rapidity gap that grows like ln(l/xp). The power law associated 
with pomeron exchange gives a distribution of events that is approximately (ixp/xp, and 
hence approximately a constant distribution in the size of the rapidity gap. That is, it gives 
non-exponentially suppressed rapidity gaps. 

In contrast, if one has a statistical distribution of particles, uniform in rapidity, and if all 
the correlations between particles are local in rapidity, then the probability of a rapidity gap 
of size Ay is proportional to e~"^^, for some constant a. Experimental data, for example. 
Fig. 1^ indicate that for small gaps there is indeed an approximately exponential fall-off with 
increasing Ay. But for large enough gaps, more than 2 or 3 units, a presumably different 
mechanism of particle production becomes significant, and there is no longer a strong expo- 
nential drop. Of course, once Ay gets close to the kinematic limit, the distribution in Ay 
must decrease again. 

Strictly speaking, one should not talk about specific events as being diffractive or non- 
diffractive. Instead one should talk about the diffractive and non- diffractive components of 
a cross-section. When one says that a specific event is diffractive, one really means that it is 
in a kinematic region (of xp or rapidity gap) where diffraction is the dominant contribution 
to the cross section. 
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FIG. 2. Distribution of 

??max ill photoproduction at HI There are no detected particles with 
pseudo-rapidity between r/max at the end of the detector at rj = 3.65 

B. Existence of a rapidity gap 

Now the maximum value of the + component of the momentum of a hadron in the X 
part of the final state in (|16]) is A"*" = Xpp'^. It then follows from Eq. that the maximum 
possible value of rapidity for any particle in X is 

2xpE 

1/kin. max ; 

m 

where E is the beam energy (assumed to be much larger than the proton mass). 

For a numerical example, consider a diffractive process with an incoming proton of energy 
900 GeV (as at the Fermilab Tevatron), and with xp = 0.01, a value typically considered to 
be safely in the diffractive region. The incoming proton itself has rapidity 7.6. From Eq. 
(^) the maximum possible rapidity of pions in the event is 4.9. There is therefore a rapidity 
gap on the proton side of the event. Since there are normally many particles in the event, 
which share the momentum, the edge of the rapidity gap, ?/max, will typically be a unit or 
more lower. 

The distribution of the soft particles in a hadronic collision is approximately uniform 
in rapidity. So when the final-state particles are plotted in azimuth and rapidity y, a 
diffractive event therefore looks like Fig. ^. When xp is decreased, the value of ?/kin. max also 
decreases, and hence the size of the rapidity gap increases. 

C. Computation of xp from the non-diffracted part of the final state 

In principle, xp can be either measured from the final state hadron: xp = 1 — p'^ /p'^ or 
from the rest of the final state: 

7 




(20) 







a' 



jet1 



-t4 



Jet 2 




2r 



FIG. 3. Plot of the distribution of final-state particles in a typical diffractive event is uniform 
in (j) and y except for an isolated forward proton. 
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(21) 



where the sum is over all particles in the remaining part X of the final state. 

To measure xp from the diffracted hadron requires first that there be elements of the 
detector (e.g., Roman pots) close to the beam-pipe, which is not always the case. Moreover, 
if is very small, the hadron's energy needs to be known very precisely to obtain a useful 
accuracy on x^; if the measurement is not sufficiently accurate, one cannot distinguish the 
value from zero. (Values of xp in the range 10"'* to 10~^ are often discussed.) So if a useful 
measurement of xp is not available, one must measure the rest of the hadronic final-state, 
and deduce xp by using Eq. (|2l|). 

In principle, the use of Eq. ( pl| ) requires knowledge of all the (typically many) particles in 
X. One might imagine that poor information on X, for example from a lack of calorimetry in 
the backward part of the detector, would prevent such a measurement from being practical. 

In fact, as I will now show, the value of the right-hand-side of Eq. (^1]) is dominated by 
the particles of the largest kr (i.e., jet particles), and by the particles that are closest in 
rapidity to the diffracted hadron, i.e., closest to the edge of the rapidity gap. 

We simply write each in Eq. ( pi] ) in terms of rapidity and transverse momentum to 
obtain 

ieX 

Here I have approximated by E\/2, with E being the energy of the incoming proton, 
which is assumed to be large. The transverse momentum appears in form of the transverse 



energy = y kfj, + . We next recall that particles in a hadronic scattering event are 
typically uniformly distributed in rapidity between some kinematic limits, and have a steeply 
falling spectrum in transverse momentum. The exponential of rapidity in Eq. (|22|) therefore 
implies that the sum is dominated by those particles of the highest rapidity. 
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The exception to this statement is given by particles that result from a hard scattering 
and that therefore have large transverse momentum. Hence the general statement is that 
the sum in Eq. ( |22| ) is dominated by (a) particles with the largest rapidity (i.e., closest to 
the edge of the rapidity gap), and (b) particles with high transverse momentum. This means 
that the measurement can be effectively be made from 

^F- (23) 

particles with highest y 
particles with large Eijp 

Hence a lack of detector coverage on the opposite side to the diffracted proton does not 
greatly affect the accuracy of the measurement. 

Exactly this method is used in a corresponding problem in deep-inelastic scattering, 
where it is known as the Jacquet-Blondel method. There it is desired to measure the 
scaling variable y (which has no relation to rapidity). This variable is defined in terms of 
the scattered lepton in exactly the way in which our is defined in terms of the diffracted 
proton. By momentum conservation, y can be computed from the hadronic final state. 
The Jacquet-Blondel method is in regular use in deep-inelastic experiments at the HERA 
collider, and it provides useful accuracy. It is therefore evident that the same method applied 
to diffraction may well provide a practical method for measuring xp.0 



D. Scaling variable /3h relative to pomeron 

When one has not only diffraction but also a hard scattering, there is another important 
kinematic variable that can be measured by a similar method. A closely related variable 
is the one called (3 in diffractive deep-inelastic scattering. I have added a subscript 'H' to 
distinguish my definition, with 'H' denoted 'hard scattering'. The definition and its physical 
interpretation can be explained by considering a diffractive hard collision treated by the 
QCD factorization formula, as formulated by Ingelman and Schlein The situation is 
represented in Fig. ^. There the hard part (i.e., the large pt subprocess) of the collision is 
treated as being due to a collision of a parton from the antiproton with a parton out of the 
exchanged pomeron. The parton in the pomeron has a distribution that is a function of its 
momentum fraction, /3h, in the pomeron. Now fractions here are defined to be fractions of 
the -|- component of momenta. So we can measure /3h from the hadronic final state by 

/Sh = ^, (24) 

Xp 

where xjet is defined in the same way as xp, except that the sum is restricted to the particles 
that are generated by the hard scattering subprocess, i.e., the large transverse momentum 
particles. Thus 



^ One counterargument is that when one measures y in deep-inelastic scattering, the particles 
dominating the sum corresponding to the one in Eq. ( p3|) are those with large transverse momen- 
tum, whereas in diffractive hadron scattering, the soft particles, those in the underlying event, 
are important. In diffractive hadron scattering one will measure all the relevant low transverse 
momentum particles only if the tracking is very good. 
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The sum in the denominator is effectively over just the particles with the largest rapidity 
and/or transverse energy. The importance of measuring cross sections differential in /?h is 
that it gives a direct probe of the parton densities in the pomeron. 

I use the name 'jet' in the symbol Xjet, in view of the fact that an important application is 
to events containing high jets. In that case the sums are actually over hadrons. However, 
it is important to observe that the same principles apply to any other hard process, like 
Drell-Yan. Then the sums are over all particles of the specified kinematics, hadronic or 
leptonic. 

Important note about notation: In diffractive deep-inelastic scattering, it is common to 
define a variable /? = Xbj/xp. In an event that is generated by the lowest-order mechanism 
of elastic electron- quark scattering, it can readily be verified that /3 = /3h. This is in fact a 
motivation for the definition of /?. 

Moreover, although the motivation for the definition of /3h, and many of its uses, are 
associated with the factorization formula (the 'Ingelman-Schlein model'), the reader should 
observe that the definition of /5h is completely independent of the model. The expression for 
/3h involves only properties of measurable hadrons (and leptons) in the collision. It does not 
involve the unobserved partons that are the mechanism by which the process is supposed 
to be generated. This observation is particularly important in the light of the fact that the 
factorization formula for diffractive hard scattering has not been proved from QCD, unlike 
the related formula for inclusive hard scattering. 



E. Relation between rapidity gaps and xp 

Much diffractive data is presented as a function of the position of the edge of the pseudo- 
rapidity gap rather than as a function of xp. So it is useful to obtain an approximate relation 
between the two variables, in order to be able to compare different sets of data. I will present 
results both for the gaps in rapidity and in pseudo-rapidity. 

1. Estimation of the size of the gap 

In a typical event, the distribution of particles as a function of y is supposed to look like 
Fig. 11(a). To estimate the relation between the edge of the rapidity gap and xp, I approx- 
imate the true distribution by a theta function distribution. Fig. ^(b). In the approximate 
distribution we have a uniform distribution, for values of y less than a value I will call yedge, 
and no particles above yedge- The height of the distribution I will call dN/dy. From the 
distribution of hadrons in real diffractive events 0]0 1 estimate that an appropriate value for 



^ There is also earlier data from the ISR H. 
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FIG. 4. Diffractive hard scattering. A parton inside the pomeron undergoes a scattering scat- 
tering with a parton out of the antiproton. 



yedge is about one unit less than the actual end of the rapidity gap, in order to obtain the 

same value for xp. 

Next, 1 split the computation of xp into a sum over the highp^ particles (e.g., associated 
with jets), and a sum over the low px particles (which comprise what is often termed the 
'underlying event'). For the particles in the underlying event, I replace their transverse 
energy by its mean value {Et). This gives 



~ (5nXf + dye 

J —oo 



dy 2E 

dN (Et) 



dy 2E 



(26) 



Here, I have extended the integral over y to — oo, since, in the underlying event, only the 
hadrons closest to the edge of the rapidity gap are important in the calculation. Hence 



In — ~ — In 



1 ^ dN 2E 

Z -T- - Vedge - In — + In — — 

1 - Ph dy {Et) 



-In 



1-/9h 



+ 



, 2E , dN 



(27) 



with the term in square brackets being approximately independent of the event's kinematic 
parameters, ymax and Pu- 



ll 



^max 

(a) 



(b) 

FIG. 5. (a) True distribution of particles in a diffractive event, (b) Theta-function approxima- 
tion to this distribution. 

2. Numerical values 

For a relation between Xf> and the maximum measured j)setidorapidity, let us rewrite Eq. 

m as 

In— = C- w-ln-^+ln— (28) 

1 — Ph 900 GeV 

with 

, 1800 GeV , dN 
C = rjmi,^ - yedge + In — 7-=r- In — . (29) 

{Et} dy 

Now C is a pure number that should be approximately independent of any of the parameters 
of the collision. So we can try to estimate it from the CDF data on single diffraction [Q. 
Their plots for the number distribution of charged tracks as a function of pseudo-rapidity 
are averaged over all events. Since different events have different values of xp, the plot need 
not exhibit a rapidity plateau even if individual events do so. Moreover the rise of dN/dy 
from zero to its plateau will be spread out. 

From the data, I estimate that at i/i = 1800 GeV 

~ 2 

drj ' 

'7max 2/edge ~ 1- (^O) 
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In addition I use the value {Et) ~ 0.3 GeV. These values give C = 8.5. 

Eq. (^81) is, of course, not exact on an event-by-event basis. However, it should be useful 
on average, with the real data being smeared around the values related by the formula, and 
of course the numerical value '8.5' should be adjusted by a fit to ordinary single diffraction 
data. 

There is indeed such data. CDF in Ref. give a plot of T^max against xp for their 
diffractive dijet data. Unfortunately, it disagrees quite badly with the above formula, both 
as to normalization and slope. This needs to be investigated. For example, the numbers I 
quoted above in Eq. (RO) may be wrong, and my formulae have made no allowance for the 
fact that the measured rj^ax applies to particles above some value of Ex ~ 400 MeV. The 
true ?7niax may be rather larger. 
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